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Abstract: We compute the spectra of the Hessian matrix, Pi, and the matrix 
M. that governs the critical point eqnation of the low-energy effective snpergravity, 
as a fnnction of the complex structure and axio-dilaton moduli space in type IIB 
flux compactifications at large complex structure. We find both spectra analytically 
in an hj -|- 3 real-dimensional subspace of the moduli space, and show that they 
exhibit a universal structure with highly degenerate eigenvalues, independently of 
the choice of flux, the details of the compactification geometry, and the number of 
complex structure moduli. In this subspace, the spectrum of the Hessian matrix 
contains no tachyons, but there are also no critical points. We show numerically 
that the spectra of PL and M. remain highly peaked over a large fraction of the 
sampled moduli space of explicit Calabi-Yau compactihcations with 2 to 5 complex 
structure moduli. In these models, the scale of the supersymmetric contribution to 
the scalar masses is strongly linearly correlated with the value of the superpotential 
over almost the entire moduli space, with particularly strong correlations arising 
for < 1. We contrast these results with the expectations from the much-used 
continuous flux approximation, and comment on the applicability of Random Matrix 
Theory to the statistical modelling of the string theory landscape. 
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1 Introduction 

Flux compactiflcations provide a promising framework for connecting string theory 
with the models and phenomena of particle physics and cosmology [1-6] . Fluxes give 
rise to an energy density that depends on the shape, and hence the moduli, of the 
compactiflcation manifold, and the minimisation of this energy allows for a controlled 
tree-level stabilisation of several moduli, which is helpful in fixing moduli-dependent 
coupling constants in the low-energy theory and in bringing it into agreement with 
observational constraints on massless scalars. Furthermore, the fluxes backreact 
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on the geometry and may create regions of the extra dimensions with signihcant 
warping. As there is a very large number of possible compactihcation manifolds and 
flux choices, the number of low-energy, four-dimensional effective theories arising 
from flux compactifications is enormous and prompts the notion of a ‘flux landscape’ 
of effective theories (see e.g. [7-11] for reviews). 

The existence of such a landscape raises several challenges. First, deriving the 
explicit predictions from generic flux compactifications with tens or hundreds of 
moduli helds is computationally prohibitively complicated: the explicit form of the 
(classical) flux-induced effective potential has so far only been found for examples 
with a handful of moduli (for which the period vector has been explicitly computed), 
and a systematic classification of the vacua for a given compactihcation manifold has 
only been possible for comparatively simple examples with few moduli (see e.g. [12, 
13]). With generic vacua out of reach, the generally applicable lessons from these 
special, explicit constructions may at best be inferred by extrapolation, which may 
be tenuous. 

Second, it is a priori possible that the landscape of hux vacua contains a large 
number of solutions that are compatible with the outcomes of any experiments and 
observations that humankind may ever conduct. If so, it is unlikely that explicitly 
detailing the properties of any given hux vacuum will lead to profound insights, 
even if the computational obstacles for constructing explicit, generic hux vacua were 
overcome. 

Rather, both these challenges motivate a statistical approach: by approximating 
the quantised huxes as continuous variables, a great deal has been learned about the 
distribution of hux vacua without requiring the direct construction and enumeration 
of the corresponding solutions [14-17] (see also [7, 8, 11] and references therein). 
The employment of statistical tools opens up the possibility of hnding limits of the 
theory where the relevant distributions take relatively simple forms due to some 
form of central-limit-theorem type of behaviour. The spectra of matrix ensembles 
provide a particularly compelling target in this respect since for large matrices with 
randomly distributed entries, the spectra quickly approach ‘universal’ limits that 
are largely independent of the statistical input. For example, the spectrum of large 
Hermitian matrices with statistically independent, normally distributed entries is 
famously given by the Wigner semi-circle law [18, 19], but so is the spectrum of 
any Hermitian matrix ensemble with independent and identically distributed non- 
Gaussian entries as long as the moments of the distribution are sufficiently bounded 
[20, 21], and so is the spectrum of random Hermitian matrices with a large number 
of statistically correlated matrix entries [22-24]. The existence of strong universality 
theorems has motivated the employment of random matrix theory (RMT) techniques 
in the study of the flux landscape [8, 25-28] (see also [29-37] for additional string 
theory motivated applications of RMT to cosmology), thus potentially providing a 
significant extension of the computational reach of the statistical methods. While the 
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relevant matrices arising in flux compactifications are not generic, random matrices, 
but rattier carry a lot of structure inherited from the geometry and topology of the 
compactihcation, one may hope that for sufficiently large and complex systems, RMT 
universality will dominate the string theory correlations, and comparatively simple 
spectra may emerge as a result [8, 25]. 

Two matrix ensembles play particularly prominent roles in the study of flux 
vacua: the hrst matrix, which we denote by A^, is formed from the second covariant 
derivatives of the superpotential W as. 


M 


f_0 

0 J ’ 


( 1 . 1 ) 


where = B^D,W = - T^^F, for F, = D,W = daW + K^W 

and Ka = daK. Here i? denotes the argument of the superpotential, and we have 
set Mpi = 1. The importance of M. is two-fold: the elements of Zab set the 
scale of the supersymmetric contribution to the masses of the chiral helds, and 
A4. appears in the critical point equation, daV = 0, of the F-term scalar potential, 
V = e^ {K°^FaFi - 3|lT|2j, as [16], 


MF± = ±2|1F|F±, 


( 1 . 2 ) 


where. 


/ pa^ I'd \ 

y±papi) ) ■ 


(1.3) 


One of the main results of [16] was to note that the symmetries of Ai are just those 
of the symmetry class Cl in the classihcation of mesoscopic Hamiltonians by Altland 
and Zirnbauer [38]. It was furthermore argued in [16] that for sufficiently generic 
compactihcations with many moduli, the rough features of the spectrum should be 
well-described by the corresponding random matrix theory ensemble. 

The second matrix of considerable interest for the counting of flux vacua is the 
Hessian matrix dehned by. 


(^a^bV 

VVaVsH VaVbHy' ■ 


(1.4) 


For a critical point to be a metastable vacuum, the spectrum of "H must be positive 
dehnite. Reference [25] (by one of the present authors, and collaborators) showed 
that for a ‘random supergravity’ in which the superpotential and Kahler potential 
are random functions, Ti is well-described by a random matrix model consisting 
of the sum of a Wigner matrix with two Wishart-type matrices. The spectrum of 
this ‘WWW-model’ was obtained analytically by freely convolving the spectra of 
the independently contributing matrices (cf. equation (4.8) of [25]). In this model, 
typical critical points are unstable saddle-points and "H has a signihcant fraction 
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of negative eigenvalues. Since the assumptions of the random supergravity closely 
matched those proposed for the flux landscape in [16], one may expect these random 
matrix results to be applicable for the complex structure and axio-dilaton sector of 
the flux landscape. 

Clearly, it is important to verify the applicability of RMT techniques to the flux 
landscape, but doing so is hard for obvious reasons: the universal limits are expected 
to be applicable precisely for the large and generic systems that are the most challen¬ 
ging to construct explicitly. The purpose of this paper is to compute the spectra of 
M and ?{ in the flux landscape, focussing on the axio-dilaton and complex structure 
moduli sector of type IIB flux compactihcations in the large complex structure limit. 
Our main results are as follows: 

• We numerically show that the spectra of M and ?{ in explicit flux compac¬ 
tihcations with two to hve complex structure moduli differ signihcantly from 
the ‘universal’ RMT spectra (cf. section 3). The string theory spectra exhibit 
strong peaks that are absent in the random matrix theory models. 

• For > 0 complex structure moduli, we analytically compute the spectra 
of A4 and "H in an h_’ -|- 3 real-dimensional subspace, S, of the full moduli 
space (cf. section 4). The eigenvalues of the matrix M come in opposite sign 
pairs and, in this subspace, the positive branch of the spectrum is given by 

degenerate eigenvalues equal to |kF| and a single eigenvalue equal to 3|hF|. 
The spectrum of "H is given by zero modes, -|- 1 degenerate eigenvalues 
equal to 27713/2 a single eigenvalue equal to 87773 / 2 . Here, as in the rest of 
this paper, 7773/2 denotes the gravitino mass restricted to the complex structure 
and axio-dilaton sector. These results hold for all hux compactihcations in 
the large complex structure limit, i.e. independently of the number of complex 
structure moduli, the ‘Yukawa couplings’ not all vanishing), and the (not 
all vanishing) hux conhguration. In this sense these string theory spectra are 
universal, albeit very diherent from the random matrix theory expectations. 
See Figure 1 for a comparison of the RMT spectra with the analytic string 
theory spectra in S. 

• We show that there are no critical points - neither supersymmetric nor non- 

supersymmetric - in S. We furthermore show that the slow-roll parameters 
are universally given by e = 4, = 8 (see section 5.1). 

• We show numerically that the scale of the supersymmetric masses exhibits a 
strong positive linear correlation with the value of the superpotential in the 
explicit compactihcations we consider (cf. section 5.3). The continuous hux 
approximation that underpins many statistical results on the hux landscape 
predicts a vanishing or negative correlation between these quantities within a 


4 




(a) Spectrum of A4 in units of \W\ (b) Spectrum of H in units of TO 3/2 


Figure 1. Random Matrix Theory spectra in black, string theory spectra in S in orange. 
Normalisation is arbitrary, and the relative heights of the delta-function peaks are taken 
to schematically indicate the degeneracy of each eigenvalue. 

broad set of assumptions, and we discuss the break-down of this approximation 
(cf. section 5.4). The strong correlation significantly reduces the frequency of 
compactihcations with large flux induced hierarchies in this region of the moduli 
space. 

• We consider flux compactihcations beyond the large complex structure limit 
and show that existing universality theorems of random matrix theory do not by 
default apply to these compactihcations. We suggest that RMT techniques may 
nevertheless be applicable to more general compactihcations if the geometric 
correlations that arise in string theory compactihcations are taken into account 
(cf. section 5.5). 

These hndings lead to several directions of possibly very interesting future research, 
and we briehy discuss these together with our conclusions in section 6. 

2 Type IIB flux compactiflcations 

In this section, we briehy review the structure of the four-dimensional supergravities 
that arise as the low-energy limit of type IIB string theory on Calabi-Yau orientifolds 
with non-trivial R-R and NS-NS hux. 

The relevant low-energy degrees of freedom for a compactihcation on the orienti- 
fold M 3 of the Calabi-Yau three-fold M 3 include the axio-dilaton, t = Co + the 
complex structure moduli, «*, where z = 1,..., and the Kahler moduli, T*”, 

where r = 1,..., To leading order in the Qs and a' expansions, the Kahler 

^ In addition, the spectrum may contain axion multiplets G“, with 0 = 1,..., which 

we will not consider in this paper. 
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potential is given by, 


K = — In J A — In {—i{T — r)) — 2 In V , (2-1) 

where V denotes the compactihcation volume and the holomorphic three-form. 
We will throughout this paper consider r in the fundamental region of the torus, 
{r G C, |Re(r)| < |r| > l,Im(r) > 0}. 

The complex structure moduli arise from the periods of as follows: for I = 
0,... take to be a canonical homology basis of and {ai,(3^) 

the dual cohomology basis satisfying. 


'M 3 


ai A A a/ = 5/ 


aiAaj= /3^A/3'^ = 0. (2.2) 


'M 3 


'M 3 


'Ms 


With respect to this basis, the periods of kl are given by. 


n = 


Iai ^ 



(2.3) 


Here serve as projective coordinates on the complex structure moduli space. The 
periods Qj satisfy the equation, 2Qj = di{z'^Qj), so that Qj is the gradient of a 
homogeneous function of degree two: Qj = djQ. 

The complex structure moduli can be expressed as the inhomogeneous coordin¬ 
ates on the space of complex structure deformations, = z^/z^ for i = 1 ,..., h^’^. 
Upon setting 2 ;° = 1, the period vector is given by. 


( 


n = 


1 

2F - v?F, 


\ 


\ 


(2.4) 


/ 


with the prepotential F = (z°) "^Q. We will be particularly interested in the ‘large 
complex structure expansion’, in which F is given by. 


F = - - KijF'U^ + XjM* -I- Kq + I . 

6 2 2 


(2.5) 


where I denotes quantum instanton contributions (that we will be more specihc about 
when considering explicit examples). In this limit, the expansion coefficients are 
given by the classical intersections of the mirror-dual Calabi-Yau, and the coefficients 
Xijk = Jms^ traditionally referred to as the ‘Yukawa couplings’. The 

d = 4, Af = 1 low-energy supergravity is obtained by orientifolding M 3 to M 3 and 
leaves the involution-odd complex structure moduli in the chiral spectrum [39] . The 
details of this involution are not important to our general discussion, and we will 
henceforth take i to run from 1 to 
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The complex structure dependent part of the Kahler potential may now be writ¬ 
ten as, 


Kr 


- In y A = - In s n) 

- In - u"){u^ - u^){u’^ - #) - 2Im(xo)^ , 


where S denotes the symplectic matrix, 



( 2 . 6 ) 


(2,7) 


We will denote the joint moduli space of the axio-dilaton and complex structure 
moduli by C. 

We are interested in compactihcations in which integrally quantised RR (T 3 ) and 
NS-NS (hfa) fluxes wrap some non-trivial three-cycles of M 3 , 


1 

(27r)^Q;' 





(27r)^Q;' 


M 3 




It is convenient to introduce the complex three-form flux G 3 
the complexihed flux vector (without subscript) as. 


As-ns e , 

( 2 , 8 ) 

= F 3 — rMs, and dehne 


TV = S 



(2,9) 


The fluxes contribute to the D3-charge tadpole by. 


Qflux 


(27r)4(a')2 


M 3 A M 3 — 


'M 3 


(27r)"‘(Q;')^ r — r 


N^EN. 


( 2 . 10 ) 


Requiring that the total sum of D3 charge vanish in the internal space leads to a 
joint condition on the D3-brane content, the fluxes, and the D7-brane and 0-plane 
conhguration, 

Qflux + Md 3 = ^ —-iVo 3 , ( 2 - 11 ) 

where iVD 3 , Nq^ denote the net number of D3-branes and 03-planes, and y is the 
Euler characteristic of the F-theory four-fold that corresponds to the given D7-brane 
and 07-plane content. Fluxes that preserve some supersymmetry contribute posit¬ 
ively to the tadpole, and when considering specihc examples we will impose. 


0 M Qfiux M ) 


( 2 . 12 ) 


where L* denotes the model-dependent maximal contribution of the fluxes to the 
D3-tadpole, given a certain conhgnration of D7-branes and orientifold planes. The 
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tadpole condition does not ensnre snpersymmetry, and flnx satisfying (2.12) is gen- 
erically non-snpersymmetric. 

The flnxes induce a complex structure and axio-dilaton dependent energy density 
that in the four-dimensional theory is captured by the flux induced superpotential 

[40], 

W = [ G3An = N-fi. (2.13) 

JM3 

Clearly, W is linear in r and, classically, at most cubic in the complex structure 
moduli in the large complex structure expansion. 

We may now dehne the ‘flux landscape’ as the ensemble of four-dimensional 
M = 1 supergravities with a Kahler potential of the form (2.1) and a superpotential 
of the form (2.13) for the set of viable compactihcation manifolds and consistent 
choices of flux. 

Note that we do not restrict the flux to be supersymmetric (as we want to study 
the moduli space dependence of the spectra of M. and H), so that the tadpole con¬ 
dition (2.12) does not bound the flux choices to be hnite for a given compactihcation 
manifold and hux tadpole. However, we will in addition require that the Kahler 
sector is stabilised in such a way that the four-dimensional supergravity is consist¬ 
ently the controlled low-energy limit of the corresponding string compactihcation. 
In other words, while we will not explicitly consider Kahler moduli stabilisation in 
this paper, we will assume that the compactihcation volume is stabilised at a suffi¬ 
ciently large value to justify the a'-expansion, the neglect of higher KK-modes, and 
the validity of the supergravity action at energies below the KK-scale. For example, 
we implicitly require that rn 3 / 2 \iu\\/'mKK -C 1, where m 3 / 2 |fuii denotes the gravitino 
mass including the Kahler moduli. Phrased in terms of the gravitino mass of the 
truncated axio-dilaton and complex structure sector, 

m3/2 = , (2.14) 

the corresponding constraint is given by [41] 

m3/2<Vl/^ (2.15) 

This ehectively bounds the ensemble of hux choices for any given manifold to be 
hnite. 

We close this section by making one additional comment on our neglect of Kahler 
moduli, as we will in this paper predominantly consider the truncated system of the 
axio-dilaton and the complex structure moduli. This truncation is well-motivated as 
we aim to study the intrinsic structure and randomness of the hux superpotential, 
which does not depend on the Kahler moduli. While the inclusion of Kahler moduli 
may change the spectra of Ti and AA through non-vanishing cross-terms, such changes 
are in many interesting cases small or easy to take into account, such as e.g. in the case 


of complex structure moduli and the axio-dilaton being stabilised at a hierarchically 
higher scale than the Kahler moduli [5], or when the no-scale symmetry of the Kahler 
sector is only weakly broken [42-46]. Moreover, not very much is known about 
the distribution of non-perturbative effects that may stabilise the Kahler moduli, 
and statistical modelling based on e.g. random matrix theory is less well-motivated 
(see however [31, 33, 35] for some interesting developments in this direction). For 
similar reasons, several previous statistical studies of flux vacua have neglected Kahler 
moduli [14-16]. 

3 Spectra of A4 and 9{ in an explicit flux compactiflcation 

We begin our study of the spectra of At and Ti with an instructive example of flux 
compactihcations on a particular orientifold. We will hnd it useful to return to this 
example at several points throughout this paper, and we will refer to it as ‘Model 
1’. In this section, we numerically compute the spectra of At and "H for canonically 
normalised helds in Model 1 as a function of the effectively four complex-dimensional 
moduli space C. 

As any given example manifold may have particularities that could bias the 
results, we will in Appendix A compute the spectra in four additional flux compac¬ 
tihcations as well as in a non-trivial modihcation of Model 1. All compactihcation 
manifolds that we consider in this paper can be constructed using toric geometry, 
and all have been previously studied in the literature: our Models 1-4 are taken 
directly from reference [47], and our Model 5 is given by the degree 18 hypersurface 
in j^ g g, which has previously been studied in [13, 42, 48, 49]. 

3.1 An explicit Calabi-Yau orientifold compactification 

Model 1 is constructed through compactihcation on one of a mirror-dual pair of 
Calabi-Yau hypersurfaces in certain four-dimensional toric varieties [47]. The com¬ 
pactihcation admits a discrete F = Zg action acting on the periods, and by con¬ 
sidering only the subsector that is invariant under this action, the corresponding 
low-energy ehective theory includes four complex structure moduli (non-invariant 
moduli can be shown to be hxed supersymmetrically). The classical prepotential is 
obtained from the intersection numbers of the mirror pair, and is in the large complex 
structure limit given by. 


Fcl. = +3 UiU4 + |m2M4 + |m 3«4 + ^^4 + +^2+^3 + f|«4 “ <(3) 

— |uiM4 — 3 U 1 U 2 U 4 — 3 U 1 U 3 U 4 — 3 U 2 U 3 U 4 — |uiM4 — 3u2uI — 3u^u\ — |u| , (3.1) 

where ^ denotes the Riemann zeta function. World-sheet instantons correct the 
prepotential at the non-perturbative level, and the leading contributions are given 
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Figure 2. Empirical eigenvalue densities of M in units of \ W\. 


by [47], 

J — 2g2i7ni3 _|_ ]^^^g4i7™4 

+ + (3 2) 

With the particular D7-brane and 0-plane conhguration considered in [47], U 2 and 
M 3 are related by the orientifold involution and only fluxes that are symmetric under 
U 2 y-)- M 3 may consistently be turned on. Thus, this model has M 2 = M 3 and effectively 
three complex structure moduli. The flux tadpole of equation (2.12) is given by 
L* = 22. 

3.2 The spectrum of M. 

We are now interested in characterising the spectra of Af and "H as a function of 
the moduli space C. Throughout this paper, we consider the spectra of canonically 
normalised fields, and for our numerical study of Model 1, we include the world-sheet 
instantons up to the second order, cf. equation (3.2). 

We begin by computing the spectrum of at a hxed point in the moduli space, 
r = M* = 10* for i = 1,..., 4, while scanning over flux integers randomly chosen in 
the range [—5, 5] with a uniform distribution, subject to the tadpole condition (2.12). 
For reasons that will become clear in section 4, we exclude the case when both the 
RR and NS-NS flux simultaneously vanish on the cycle whose period has cubic terms 
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in the complex structure moduli (cf. the cycle appearing in (2.3)). The resulting 
eigenvalue histogram is shown in Figure 2(a), when plotted in units of |hF| for each 
flux choice. The spectrum shows prominent peaks at ±|hF|, and smaller peaks at 
±3|kF|, in stark contrast with the smooth ‘semi-circle-like’ Altland-Zirnbauer Cl 
(AZ-CI) spectrum of Figure 1(a). 

To investigate the dependence of the spectrum on the values of the moduli helds, 
we have performed various joint scans of subspaces of the moduli space and of the 
flux numbers. In Figures 2(b)-2(d) we plot the resulting eigenvalue densities for three 
cases in which complex structure moduli vevs are (for simplicity) sampled uniformly 
in the range, 

Re(n*) e [-10,10], Im(n*) e [1, , (3-3) 

with G {10,5,2}. The axio-dilaton is simultaneously sampled uniformly in the 
fundamental domain with 2 <Im(r)< where G {10,5,5}. 

The shape of this spectrum and the implications that follow from it are the main 
themes of this paper. Here, we will merely make the following simple observations: 
i) none of the densities resemble the AZ-CI spectrum, ii) the peak at ±|IF| is visible 
in all cases, while the peak at ±3|fF| is clearly visible in all but the last case. The 
peaks are less blurred for larger typical values of the moduli, in) The peaks are only 
visible in the spectrum when plotted in units of |IF| (which varies from realisation 
to realisation) and for canonically normalised helds. This could possibly explain why 
this effect has not been previously observed in the literature. 

3.3 The spectrum of "H 

The eigenvalues of the Hessian matrix for canonically normalised helds give the 
squared physical masses of the scalar helds in the theory. We now compute the 
spectrum of l-i in Model 1 in the same cases considered above for the spectrum 
of M.. Figure 3(a) shows the eigenvalue density of "H when scanning over huxes 
at r = n* = lOh In stark contrast to the random matrix theory spectrum of the 
‘WWW’ model of Figure 1(b), the spectrum exhibits sharp peaks at 0 and 
and a smaller bump at 

We furthermore consider joint scans of huxes and moduli vevs within the sub¬ 
spaces dehned as per the discussion around equation (3.3). The resulting spectra are 
shown in hgures 3(b)-3(d). The peaks in the spectra remain distinctive and promin¬ 
ent in the hrst of these cases, while they blur as the sampling is restricted to regions 
with smaller moduli values. In the hnal case plotted in Figure 3(d), the peaks have 
been blurred into a broad feature that peaks around zero. 

We make the following observations: i) In none of the cases is the spectrum 
well-described by the ‘WWW’ random matrix model of Figure 1(b). ii) The peaks 
in the spectra are somewhat less sharp than those observed in the spectra of Ai for 
the same regions, but are similarly sharpened for large moduli values. Hi) Again 
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Figure 3. Empirical eigenvalue densities of T-L in units of 

we note that the peaks would not appear very prominently had we not canonically 
normalised the helds and expressed the histogram in units of ^ 3 ^ 2 - 

In sum, we have in this section shown that in a particular, explicit flux com- 
pactihcation, the spectra of M. and 'K show curious peaks that become particularly 
prominent at large complex structure, but that appear to influence the spectrum 
over the entire large complex structure expansion. For both M. and H, we found the 
spectra to differ significantly from those of the corresponding random matrix theory 
models. We will now argue that this peaked structure is to be expected as a general 
feature of flux compactifications close to the large complex structure point, and that 
this structure has signihcant implications for moduli stabilisation and the statistics 
of flux compactifications in this region of moduli space. 

4 Analytic derivation of the spectra of A4 and Ti 

We are now interested in explaining the observed structure in the spectra analyt¬ 
ically for generic flux compactifications at large complex structure. Unfortunately, 
in general the superpotential (2.13) is quite complicated and we know of no way of 
directly diagonalising M. and H. We will however find in sections 4.2 and 4.3 that 
it is possible to compute the spectra of these matrices algebraically for general flux 
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compactifications - but only in a particular limit of the complex structure moduli 
space. We now define this subspace and briefly discuss its properties. 

4.1 The subspace S 

We will consider compactihcations with non-zero R-R or NS-NS flux along at least 
one of the three-cycles of the compactihcation manifold. Given such a compactihca- 
tion, we may choose the homogeneous coordinate to be the period of G correspond¬ 
ing to a cycle, A^, with non-vanishing flux, and we may express the inhomogeneous 
coordinates as discussed in section 2 as = u^/u^. For future convenience, we will 
denote the flux along this cycle by 7 ^ 0, without the vector arrow. In the large 
complex structure expansion, the superpotential (2.13) is then a polynomial of de¬ 
gree three in the complex structure moduli, and the cubic terms are all proportional 
to the single complex flux N. 

For sufficiently large complex structure moduli vevs, the superpotential is well- 
approximated as a homogeneous function of degree three in the complex structure 
moduli, 

W = ^Xijku'u^u^ -f 0{{u"Y) K. yXjjfcuVn''. (4.1) 

We here assume that IF 7 ^ 0 so that, in particular, the ‘Yukawa couplings’ >Cijk do not 
all vanish. We note that the exact region for which (4.1) will be a good approximation 
will depend on the fluxes on all cycles of the compactihcation manifold. The complex 
structure Kahler potential simplihes similarly at large values of the complex structure 
moduli to, 

Rc.s. ~-hi -«)*(«-m)^(m - h)^^ . (4.2) 

We will hnd it hard to make analytical progress for arbitrary phases of the 
complex structure moduli. However, the computation simplihes signihcantly in the 
subspace in which the complex structure moduli all have the same phase, i.e. 

M* = MS* , (4.3) 

1 2 

for the hj real parameters s*. We do not impose any restrictions on the value of the 
axio-dilaton, so the real dimensionality of the subspace in satisfying (4.3) is -|- 3. 
We will denote the subspace of the form (4.3) in which (4.1) and (4.2) provide good 
approximations by S. For clarity, we will label equations that only apply in 5 by a 
[iS] in the margin. 

We now introduce a short-hand notation that we hnd useful in the explicit com¬ 
putations that will follow. For some tensor Rp-fc/mn we denote tensor contractions 
with M* by subscript w, and contractions with w* — m* by A so that for example, 

Aijklmn'U {u 'U)'^'U {u m) ^uAMAmn • (^•^) 
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In this notation the superpotential and Kahler potential at large complex structure 
are given by, 


N 

W = —X, 


K = - Inl -Xaaa ) - ln(-i(r - r)) 


(4.5) 


Using (4.2), the components of the Kahler metric, = d^jK, and its inverse 


Rab given by. 


= 


r - T 




= 6 - 




IjA 


K - = K- = 0 


, ^iAA^jAA 


X, 


AAA 


X; 


AAA 


K"" = -(r-f)% ir” = ir" = 0. 


= 


-1 


ijA -{u-uy{u-uy. 


(4.6) 


We have here assumed that Xjj/\ is invertible with the inverse (in our numerical 
studies we simply ensure that the Kahler metric is positive dehnite). 

In this notation, the F-terms of the complex structure moduli are given by. 


17 n TT7 ^ ( ^UUU 

Fi I ^iuu ^i/\/\ 

4 V ^AAA 


(4.7) 


In the subspace S, the F-terms are given by, 

u 


N 

~2 


^ ^iss^ ( f 


U — U 


= -3IF 




X., 


u/ u — u 


(4,8) [5] 


where Xj^s = XjjfcS-^s^ and x^^^ = ^Cij^s^s^s^ 7 ^ 0. Using the inverse metric of equation 
(4.6), we hnd that 

F* = -Ws\u - u) Q) , (4.9) [5] 

so that FjF* = 3|IUp. Since furthermore F^- = (|)^ fU so that F^-F'^ = |IUp, the 

total contribution from the F-terms to the scalar potential is given by. 


FaF^ = F^F^ + FiF^ = A\wy . (4.10) [5] 

Thus, the moduli F-terms source a positive vacuum energy, 

V = e'f|H'f = (4.11) [SJ 

where we (still) have ignored Kahler moduli. 
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4.2 The spectrum of J\4 


We are now interested in computing the spectrum of A4 for the canonically norm¬ 
alised complex structure and axio-dilaton helds. The symmetries of M dictate that 
its eigenvalues come in opposite sign pairs, and that the positive branch is given by 
the square root of the eigenvalues of {ZZ)ab = ZacZ\. The Kahler potential (4.2) 
leads to non-trivial kinetic terms for the scalar helds at a generic point, however, we 
can obtain the spectrum of {ZZ)ab for canonically normalised helds by computing 
the spectrum of the non-Hermitian matrix K°'°‘{ZZ)ab = {ZZY^. Thus, to hnd the 
spectrum of Ai for canonically normalised helds, we compute the square-root of the 
spectrum of {ZZ)°‘i^. 

Before we go into the details of the computation of the spectrum of A4, we 
note that in hux compactihcations, the tensor Zab inherits much structure from the 
underlying Calabi-Yau geometry [15]. In particular, Zab has only independent 
complex entries, while a generic complex symmetric tensor of the same dimensions 
would have {h]f + 2){h]f + l)/2 independent components. The correlations that 
limit the number of independent degrees of freedom in Zab in string theory can in 
part be traced to the simplicity of the axio-dilaton dependence of K and W from 
which it follows that, 

Zrr = ^. (4.12) 

More signihcant however, are the correlations in the complex structure sector that 
arise directly from the Hodge decomposition of covariant derivatives of the holo- 
morphic three-form. While of course is a (3,0)-form, DiVt is (2,1) and forms a 
symplectic basis of One can furthermore show that DiDjQ is (1,2), and 

may then be expanded in terms of the basis vectors DiQ [50]. The expansion coeffi¬ 
cients are simply proportional to the ‘Yukawa couplings’ as, 

DiDjVl = . (4.13) 

We now note that the ij components of Zab can be simplihed as, 

% = DiDjW = J Gs A D.DjQ = . (4.14) 

Furthermore, the ‘mixed’ components between the axio-dilaton and the complex 
structure moduli are given by, 

z^i = Dr{N ■ An) = K^N* ■ An. (4.15) 


We thus have. 




(4.16) 
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and we may take Z^i to be the independent complex components of Zah in flux 
compactifications. The relation (4.16) is qnite nsefnl in simplifying the compntations 
in this section. 

We are now ready to compnte the spectrnm of Using eqnation (4.5) we 

And that, 


■^ri — 








K, 


iAA' 


X, 


AAA 


= -- 
^ 2 L 


^^ijm^mkA^^uu 


X, 


ijA 


X^AAA 


‘^^uuu 4“ 3x, 


Auu 


(4.17) 


We now specialise to the snbspace S and raise one index of Zab to And Z\. The rr 
component is trivial, Z'^^ = 0, and the other components are given by. 


Z^ = Wp{ 5^ 




C c o 

Z" = -3-lU^ 
■' P ^sss 

Z™ = -W—s ^, 
pc 


where we have introdnced. 


r 


c = 


u — u 

The matrix (ZZ)“^ is then given by. 


p = 


u 


2 

1 ) 


_ ^ 
^~W 


(Z7Y - \ _ ITT/P 

‘ - \(ZZ)A (ZZ)- I - 1 ’^ I 


3 Qcp^q^ 
^ 5) + 

cp-^q 3 Xsss, 


(4.18) |5] 


(4.19) 


(4.20) |5] 


Using the expression for the F-term in 5, cf. eqnations (4.8) and (4.9), we And that 
this matrix can be written as. 


(4.21) [5] 

Dedncing the spectrnm is now trivial: any linearly independent set of vectors that 
are perpendicular to Fb are eigenvectors of (ZZ)“j, with the eigenvalue \W\^. The Anal 
eigenvector is given by and since according to equation (4.10), FaF°‘ = 4|lUp, 
the corresponding eigenvalue is equal to 9|hFp. 

The spectrum of in 5 for canonically normalised flelds is thus given by. 


{ZZ)\ = |lUp5^ + 2F“Ft 


Spectrum(Wl) 


|hF| multiplicity 
3|hF| multiplicity!, 
-|fU| multiplicity , 
-3|hF| multiplicity!. 


(4.22) 

[ 5 ] 
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Equation (4.22) is one of our main results, and clearly exhibits ‘clustering’ of 
the eigenvalues into delta-function peaks in the spectrum. The location of these 
peaks are exactly those where we observed the peaks in the spectrum in the explicit 
example of section 3, cf. Figure 2(a). Thus, equation (4.22) explains the presence 
of the peaks in the observed spectrum, and moreover proves that such peaks are 
universal in the large complex structure limit: our analytical derivation applies to 
any non-vanishing choice of flux, any not all vanishing ‘Yukawa couplings’ and 
any number of complex structure moduli, 


4.3 The spectrum of H 

We now show that the spectrum of the covariant Hessian, Ti, similarly takes a very 
simple form in S. The Hessian matrix is given by. 


n = 




ab 




(4.23) 


^ ( Za Zl^ FaFl ^R^l^^F^F^ ^ Uabcr ^ 

I „ c „ 17 17 D -'pc'pd I ' 

^bc — -Tbra — ^bacd^ ^ 




0 K-abl^F^ - 2\W\^^ 


(4.24) 


where = FaF°‘, Uabc = DaD^DcW is complex and symmetric, and Rbad = 
KijjdcT^^ denotes the non-trivial components of the Riemann curvature tensor on 
the held space. We will again hnd the spectrum of R for canonically normalised 
helds by computing the spectrum of R contracted with the inverse of the Kahler 
metric, schematically ^K~^R\ 


4.3.1 The diagonal blocks of R 

Raising an index of the diagonal block matrices of R gives, 

e-^K'^^Vl^V = (ZZ)\ - FbF“ - R^-^^^F^F^ + 5^ [f'^ - 2\wf^ . (4.25) 

As we have already computed (ZZ)“^ in equation (4.21), the only non-trivial term 
remaining is the curvature term, K°‘°'R^^^cdF^■ We now show that also this contri¬ 
bution takes a very simple form in S. 

As there are no cross-terms between the axio-dilaton and the complex structure 
moduli in the Kahler potential, the curvature tensor has vanishing mixed components 
between these sectors. The axio-dilaton contribution is given by, 

K^^Iirfr,F’F’ = (T - fr (r - ff\W\^ (4.26) 

= 2|Wp. (4.27) 
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The complex structure components are not much harder to compute, and we hud 
that, 

= 6\WpSi , (4.28) [5] 

SO that, 

K^RijuF^F = K-iJ {Kyu - F^F = 2\WlFi. (4.29) [5] 

The curvature contribution to the Hessian matrix is evidently proportional to the 
unit matrix, 

= 2\W\5t. (4.30) [5] 

The diagonal blocks of K~^'H are then simply given by, 

e-KK^-a^l^Y = |pp|2^a + ^ 

e-KR-^a^l^y ^ |^|2^a ^ pap_ _ ^4 3^^ 

4.3.2 The off-diagonal blocks of 1-i 

The off-diagonal blocks of K~^'H are given by K°'^UahcF‘^ — WZ^, and its complex 
conjugate. Equation (4.18) gives the components of and we here compute the 
components of K°‘^UahcF^- The reader interested in the resulting simple expression 
- but not the intermediate technical details of the computation - may skip ahead to 
equation (4.48). 

The [/-tensor contribution to is given by, 

K^^UabcF^ = {daZbc + KaZbc - T^.Z,, - T^.Z,,) F^ . (4.32) 

For the free indices (f, r), this expression vanishes identically since D^-Zt-j. = DiZ^-^- = 


0 so that. 

K^^UrrcF^ = K^^UrrrF^ + K^^UrriF^ = 0 . 

(4.33) 

For free indices 

(f, i) we have. 



K^^UricF^ = K^^Ur^rF^ + K^^U^jiF^ = , 

(4.34) 

while for (j, r) 

we have. 



K^%rcF'' = RFUprF^ + RFUjkrF^ = R^UjkrF^ . 

(4.35) 


Clearly, the (j, r) and (r, j) components are related by, K^^UjkrF^ = [K'^'^UriiF’‘~\. 

Finally, the {j,i) components are given by, 

KFUrF^ = RFUrF^ + K^Up^F'^ . 


(4.36) 


Thus we need to compute Urij and Uijk and contract these with the relevant F-terms 
and the inverse metric. 

Using equations (4.15) and (4.13), we have that, 


Urij = DjZir = KrN* ■ DjDiYi 


In the subspace S, this expression evaluates to. 


u 


1 


Urij 6 ( _ . , -\2i -\ ’’ 

\u/ [u — uY[t — t) K. 


(4.37) 


(4.38) [5] 


Using (4.16) to simplify Ukij we have, 

Ukij = 

= {2KkK,,i% + - Vkj^^ml%). (4.39) 

We raise an index of equation (4.17) to hnd. 


z" = ^ 


^AAA l 




ikA^iuu 


uu (^ U) I Q ^uuu T „ ^Auu 


(4.40) 




To hnd dkZ we need to differentiate which is perhaps most simply done by 
noting that dk{>c'[^;,>CmnA) = 4(^0 = 0, so that dk{>t:~i\) = ->iimA^rnpk>ipnA- We 
then hnd that in 5, 


dkZ r = S^,KrW(-] (2 


M\ ' 


U/ 


(4.41) [S| 


The relevant Christohel symbols are given by, 

3 


^kj U-kjk ^IkA^^^j 


^kjm^mis + 3-5 


U — U 


^AAA 
I ^kjs 


i^i^jAA + Sj^kAA -{U- 




351^-3(55^ 


K. 


sss 


X. 


sss 


(4.42) [S| 


It is then straightforward to show that. 


KmjlZ;^ = -IF I - ] 


u\^ u — u 


uJ T — T 


-X, 


mjs ? 


(4.43) [S| 




= -W[- 


U\‘ 


X., 


-.u/ T — T 

All but one term in Ujik then cancels, and we have. 


I I c-t^jss^kis t-i^jks^iss c-,^ijs^kss \/. . r (p-i 

3 ( Xfcq + -3 ^^j3-3^^- (4.44) [S] 

^SSS J 




Ukij = -6 ( - ) 




lU 


^ijk 


uJ {u — uY X; 


(4.45) [5] 
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We are now ready to contract the components of the [/-tensor found in equations 
(4.38) and (4.45) with the F-terms. The relevant expressions are given by, 


f/ i/i = 6 (- 
\u 


1 


UijtF’’ = 61- 


(t — t)(u — u) 
2 1 

|W| 




2 ^ijs 


{u — uy 






2 1 


[u — u) 


■|W| 


2 ^ijs 


(4.46) [S| 


Upon contraction with the inverse Kahler metric we find the components, 

{K-^UF)i = 0, 


{K-^UF)j = -6 — \Wy-"" 

pq X, 

(K-^UF)l = -2 —IWI"/, 

cpq 


sss 


{K-^UPy. = 2 p\wy 51 - 3s^' 


^Xi, 


X, 


(4.47) [S| 


By comparing equation (4.47) with equation (4.18), we arrive at our final result, 

{K-^UF)l = 2WZ~^. (4.48) [5] 

Thus, the off-diagonal blocks of K~^'H are simply given by. 


-K 7 >-aav 72 


e = {K-^UF)l - WZl = WZ^ , 


(4.49) [S| 


e ■■K’^V^V = (K-^UF)f - WZf = WZf. 

4.3.3 The spectrum of "H 

Putting our results from section 4.3.1 and section 4.3.2 together, we have found that. 


K-^n = ( |fUpl + ( 


WZi 0 

= ml/2t + + e 


+ 


K 


F^Fi 

0 

0 

F^Fb 

' F^Fi 

0 

0 

F^Fb 


(4.50) [5] 


Just as for the spectrum of AJ, the spectrum of H for canonically normalised fields 
can now be read off by inspection. We hrst note that, 

70 . fpb 


Z^F^ = 3hFF“ 


(4.51) [5] 


so that the vectors F± of equation (1.3) are eigenvectors of K~^Ai with the corres¬ 
ponding eigenvalues being equal to ±3|hF|. It then immediately follows that F± are 
also eigenvectors of K~^'H with the eigenvalues. 


mo± = 


8m2/2 

2m2/2 


(4.52) 


[5] 
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The remaining eigenvectors of K are then formed from the eigenvectors of K 
with eigenvalues ±|hh|. These are all perpendicular to F± and the corresponding 
canonically normalised squared masses are given by, 



2m2/2 

0 


(4.53) 


[ 5 ] 


for i = 1,..., hh^. 

In sum, the spectrum of "H in 5 for canonically normalised helds is given by. 


(4.54) 

[ 5 ] 

Equation (4.54) together with equation (4.22) are the main analytical results in 
this paper. Just as in the spectrum of AJ, the spectrum of the Hessian matrix in 
S exhibits delta-function peaks, and we note that the location of these peaks are 
exactly as expected from the numerical analysis of Model 1, cf. Figure 3. While the 
peak at zero appears to be stronger in the numerical spectrum of Model 1, we note 
that this can be largely attributed to the different widths of the peaks around zero 
and 2m3/2- W^e thus conclude that equation (4.54) provides the analytical explanation 
of the observed peaked spectra of the Hessian matrix at large complex structure. 



[0 

1 2 

multiplicity hj , 

Spectrum('H) = < 

2m2/2 

multiplicity -t-1, 


, 8m2/2 

multiplicity 1 


5 Discussion 

The main results presented in this paper have several important - and perhaps 
unexpected - implications for scenarios of moduli stabilisation as well as for the 
statistics of flux compactihcations. In this section, we discuss what we regard as the 
most important of these implications. 

5.1 No critical points and large slow-roll parameters in S 

We begin by making a few simple observations regarding the subspace S. In section 

4.1 we found that = 4|Wp in S. Since we have assumed that \W\ ^ 0, there 
are no supersymmetric vacua in this subspace. Furthermore, the spectrum of M. in 
5, cf. equation (4.22), has no support at 2| W| so that the critical point equation (1.2) 
has no solutions. Thus, there are in S no critical points of the truncated axio-dilaton 
and complex structure moduli system. 

Furthermore, the inflationary slow-roll parameters take universal, large values in 
S. Since, 

daV = e^WF^, (5.1) [5] 
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we find that, 


1 

^ “ 2 


2daVK°^diV 


= 4, 


(6.2) [S| 


which is independent of the flnx choice and details of the compactification manifold. 
We define the slow-roll ri\\ parameter to be given by, 


^11 


V 


(5.3) 


where ca = —dAV/\\dV\\ for ||c?ld|| = ^/K^^dAVdsV, and A rnns over both holo- 
morphic and anti-holomorphic indices. Using onr expressions from section 4.3, we 
hnd that. 




= -mi,, . 


It then follows that. 


^11 


(5.4) [S| 
(6.6) [S| 


in S. We have thus shown that there are no vacua in S and that the slow-roll 
parameters are much too large to support slow-roll inflation. 

It is important to note that the inclusion of additional moduli helds may alter 
these conclusions. We illustrate this with the example of the class of approximately 
no-scale supergravities considered in [46] where it was shown that de Sitter vacua can 
be constructed in string theory motivated M = 1 supergravities through only a very 
limited amount of tuning. In this scenario, supersymmetry is predominantly broken 
by a ‘no-scale’ field (cf. Kahler modulus), with only a small amount of supersymmetry 
breaking arising from other helds (cf. the axio-dilaton and the complex structure 
moduli). The critical point equation for these ‘other’ helds can be written as an 
eigenvalue equation of the truncated matrix Wl, or perhaps simpler in terms oi ZZ 


as, 

(ZZ)/Ft=|4l'|2F., (5.6) 


for the canonically normalised helds in this sector. In other words, while the full crit¬ 
ical point equation requires M. to have an eigenvalue equal to 2|1U|, when truncated 
to the helds perpendicular to the no-scale held, M. should have an eigenvalue equal 
to |1U| to solve the critical point equation. This clearly illustrates that a study of 
the axio-dilaton and complex structure sector alone does not suffice to make general 
statements of the existence of non-supersymmetric critical points.^ Controlled mod- 

^In 5, [ZZ)^ indeed has (several) eigenvalues equal to \W\'^, however, it is not hard to see 
that the solutions of the type suggested in [46] cannot be constructed in this subspace: the F- 
terms of the complex structure moduli and the axio-dilaton are not small with respect to |1T| in 
S, and furthermore, the F-term is exactly the eigenvector of ZZ with eigenvalue equal to 9|lTp, 
cf. equation (4.51). 
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uli stabilisation is also necessary for the construction of viable models of inflation in 
string theory. For a discussion of some issues that arise in complex structure moduli 
inflation, see [51-53]. 


5.2 Highly peaked spectra at not-so-large complex structure 

We now return to our explicit example, Model 1, and ask the question: for a given 
choice of flux numbers and as a function of the moduli space, by how much does the 
spectrum deviate from the analytical spectrum that we found in section 4 for the 
subspace 5? How close does it come to the predictions of the random matrix theory 
models? 

To quantify the discrepancy between eigenvalue spectra, we introduce the ‘spec¬ 
tral deviation’, pi, as a measure of the fractional difference between two given eigen¬ 
value conhgurations. We define this measure as follows: take a = (ui,..., uat) and 
/3 = (/?i ,..., Pn) to be two sets of eigenvalues, listed in increasing order i.e. ai < 
a 2 < ■ ■ ■ < cxn, and (di < (^2 < ■ ■ ■ < I^n- The spectral deviation of (3 from a is then 
given by, 

(6.7) 

where the product is taken to be the ordinary Cartesian inner product. 

We will consider two reference vectors a when computing the spectrum of M. as 
a function of moduli space. The hrst one is given by the mean positions of the 
eigenvalues in the AZ-CI random matrix ensemble, which we hnd numerically to be 
given by, 

aAz-ci = (0.35, 0.80,1.27,1.80,2.45)|W|, (5.8) 

for five fields, where we have included only the positive branch eigenvalues. The 
corresponding AZ-CI eigenvalue density (that we will not use but include here for 
clarity and reference) is given by. 


Pci(A) 


1 

2Ti da‘^\X\ 




A2)(A2-r/_), 


(5.9) 


where ri± = da‘^{l ± a/1 + 1/d) for Ai being 2d dimensional. The reference vector 
«AZ-ci corresponds to a spectrum with cr^ = 2/5. 

The second reference vector is the positive branch of the spectrum of Ai in S, 
i.e. ds = (1,1,1,1, 3)|hF| for hve helds. To illustrate its use, we note that the spectral 
deviation of (Jaz-ci from the ‘S spectrum’ for hve helds is given by /^^(aAz-ci) = 
0.34, while PAz-ci{dls) = 0.36. 
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(c) fis- line at fis = 0.02 


01 


(d) /15: line at fis = 0.02 


Figure 4. Spectral deviations for the fluxes of equation (5.10). In Figures (a) and (b) 
= u and r = lOi. In (c) and (d), aig{ui) = 6i and |ri*| = 5 and r = lOi. In (c) 
02 = 03 = 7r/2, in (d) 04 = 7r/2. 


In Figure 4 we illustrate the moduli dependence of the spectral deviation for 
‘coincidental’ complex structure moduli, m* = m, and r = lOi and the flux choice, 


1 

{2'KYa' 

1 

(27r)^a' 



{2, 1 ,- 5 ,-5,4}, 
{5,2,5,5,2}, 


i2Tiya' 


(27r)^a' 


'Bi 



-03 — { 4, 4, 3, 3, 4} , 

= {2,-5,-1,-1,-3}, (5.10) 


that contributes to the D3-tadpole by Qhux = 3. 

In this example, the AZ-CI spectrum does not provide a good approximation 
of the spectrum in any part of the sampled moduli space, even for small imaginary 
parts of the complex structure moduli, cf. Figure 4(b). 
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The 5-spectrum is quickly approached for imaginary parts of the moduli vevs of 
C>(few), with /i 5 < 0.05 over a signihcant fraction of the moduli space, cf. Figure 4(a). 
For large complex structure moduli, the deviation is of the order of a few percent for 
general values of the moduli phases. In Figures 4(c) and 4(d) we illustrate this by 
scanning the phases of the complex structure moduli between 0.27r and O.Svr while 
keeping \ui\ = 5. 

The moduli dependence of the spectral deviation certainly depends on the fluxes, 
but the example shown in Figure 4 is not atypical. Consistent with the expectation 
from our analytical derivation, we hnd that large flux numbers on the A'^-cycle leads 
to the large regions with small /i^. In no case have we found that the AZ-CI spectrum 
provides a good approximation to the spectrum of M. over a signihcant fraction of 
the moduli space. 

5.3 Strong linear correlation between W and the snpersymmetric masses 

We now discuss one particularly interesting consequence of the strong peaks in the 
eigenvalue spectrum of M.. From section 4 we have seen that in 5, the spectrum of 
M. is given by integer multiples of \W\ so that, in particular, the spectrum of M. is 
perfectly positively linearly correlated with |1F|. On the other hand, it is frequently 
assumed in the literature that the scale of the supersymmetric masses (set by Ad), is 
statistically independent of the magnitude of the superpotential, i.e. that the above 
correlation should instead vanish. We will discuss the motivation for this assumption 
in more detail in section 5.4. In this section we numerically compute the correlation 
as a function of the moduli space for Model 1. 

The scale of the supersymmetric masses has important phenomenological con¬ 
sequences. If the supersymmetric masses for a set of moduli helds can be arranged 
to be much larger than W and the F-terms, then the Hessian matrix of this sector 
becomes positive dehnite, with the diagonal blocks given by ZZ. Such moduli are 
then supersymmetrically stabilised, and can consistently be integrated out. This 
strategy was famously employed in [5] to stabilise the complex structure moduli and 
the axio-dilaton at the flux scale which, by tuning of flux numbers, was taken to be 
signihcantly larger than the value of | 1 F|. 

The scale of M. is clearly set by the scale of Zah, but we have seen from equations 
(4.12) and (4.16) that the independent components of the Z-tensor can be taken 
to be Z^-i. The magnitudes of these components set the scale of the flux-induced 
supersymmetric masses, which we dehne as, 

ms = Zr^^" . (5.11) 

We are then interested in the correlation between and the gravitino mass, m^i 2 = 
0 -^/ 21 ppI, g^g clehned in equation (2.14). 

We emphasise that here, the prefactor exp(iF/2) is taken to include the Kahler 
potential truncated to the axio-dilaton and complex structure moduli sector, and an 
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Figure 5. Pearson’s correlation coefficient as a function of complex structure moduli with 
u* = u. Lines indicate r = {0, 0.25, 0.5, 0.75}. 

additional volume suppression of this scale will appear in the full compactihcation 
including Kahler moduli.^ 

To extract the strength of the correlation between rris and ^ 3/2 as a function of 
the moduli space, we scan over fluxes to create an ensemble of values of (rUs, 777 . 3 / 2 ) 
at each given point in the moduli space. From this ensemble, we compute Pearson’s 
correlation coefficient, 

Covariance (777^, 7773/2) 

where a denotes the relevant standard deviation. In practice, we consider ensembles 
generated from 1000 random flux choices (subject to the tadpole condition 0 < 

particular consequence of not including Kahler moduli is that TO 3/2 defined by (2.14) can 
(and typically do) take on values larger than the string scale. Such flux compactifications require 
large volume suppressions to render the EFT to be controlled, as discussed in e.g. [41]. 


(5.12) 


26 






































(a) T = 5i, r = 0.92 (b) r = lOz, r = 0.95 

Figure 6. The distribution of (m^, 7713 / 2 ) at u* = —1.6 + 1.2i for 10,000 random flux choices 
consistent with the tadpole condition. 

Qflux < 22) at each point of the densely sampled moduli space. This enables us to 
study the dependence of the correlation coefficient r on the moduli helds. In S, we 
expect perfect positive linear correlation with r = 1, while in the regions in which 
rris and 7773/2 are independent we should find r = 0. 

In Figure 5 we plot r for the ‘coincident slice’ of the complex structure moduli 
space for which = u for all i, and consider the string couplings = 1, and 
The strength of the correlation increases with decreasing string coupling, and for 
Qs < 1/5, rris and 7773/2 are strongly linearly correlated over a very large fraction of the 
sampled moduli space. We have again verihed that this result does not signihcantly 
depend on the assumption of ‘coincidence’ of the complex structure moduli vevs: 
qualitatively similar results arise even for random phases of the moduli. 

To illustrate the strength of the correlation between rris and 7773/2, in Figure 
6 we plot the distribution of pairs {rris, 1115 / 2 ) obtained from scanning over 10,000 
flux choices at a particular point in the moduli space. The strong linear correlation 
between the two quantities is plainly visible from the plot. Could this correlation be 
due to some particularity of the explicit models that we study? We will now argue 
that the answer to this question is no: a strong linear correlation should be expected 
for any compactihcation close to the large complex structure point. 

5.4 Why does the continuous flux approximation break down? 

A key tool in the derivation of many statistical results on the flux landscape - 
including the famous derivation of the index density of supersymmetric vacua in 
[14] - is the approximation of integer quantised fluxes as continuous variables. This 
approximation is often assumed to be valid for sufficiently large flux tadpoles and 
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allows for the replacement of sums over integers with continuous integrals, 


J dA^RRadiVNSa0(Qflux)0(-h* — Qflux) • (5.13) 

(.^RR . .^Ns)lo<Qfl^,^<L* 

For our purposes, it is convenient to make the following change of variables: for any 
given value of the complex structure, {Q,ri, Difl, D^Cl) form a basis of and 

we may correspondingly expand the flux vector with respect to the basis of peri¬ 
ods (n, n*, Hjll, This basis is not orthonormal in general, and the relevant 

symplectic inner products are given by, 

n s n* = , An s A-n* = . ( 5 . 14 ) 

To simplify our expressions, we phrase our results in this section in terms of the 
transformed flux vector, 

iV= j = SiV. (5.15) 

In the new basis, this transformed flux vector has the expansion, 

N = 

This is a linear change of variables from the 4(h)_’^ -|- 1) real fluxes in the canonical 
basis to the 2 -|- 2 {hh^ + 1) + real components of (IF, A, Ai)-^ The integral 
(5.13) is in this basis given by, 

2(hF+l) 

n diVRR,diV^5a = C /d^IFd^^'^-'+dir^d^'^-'A*, (5.17) 

a=l 

where C is a complex structure dependent - but IF, F, and independent - 
constant. 

Using the expansion (5.16), we see that the flux contribution to the D3-tadpole, 
cf. equation (2.10), is given by, 

[(27r)A'2] Qfl,, = e"" (jlFp + - e^AA • (5.18) 

Thus, in the continuous flux approximation the correlation between rrts and ^ 3/2 
arises solely from the tadpole condition (5.18) and not from the measure (5.17). 

^We are grateful to Kepa Sousa for discussion of this point. 

® In this basis, the familiar ‘ISD condition’ which stipulates that 3-form flux, G3, that pre¬ 
serves Af = 1 supersymmetry has vanishing (3,0) and (l,2)-components becomes quite transparent: 
supersymmetric configurations have no support along 11 and DjTi*. 


F- - Z - 

IFA - An + A'A-n* - 


(5.16) 





There are two particularly interesting cases to consider: first, for a sufficiently 
large tadpole L*, there will be numerous flux choices giving ^ 3^2 ^ [( 2 vr)^a'^]L* for 
which the tadpole condition Qaux < -h* does not introduce a correlation between 
rus and ^ 3 / 2 . In this case - which is perhaps the most commonly considered in the 
literature - we have r = 0 . 

Second, we may consider the general case in which and ^ 3/2 are random 
variables of the flux choices, and typically not much smaller than [{2'KYa']y/L^. 
We will not discuss this general case in full detail, but rather consider the simpler, 
supersymmetric case. We furthermore take 777 . 3/2 and rris to be uniformly distributed 
in the quarter disc of radius [{27i)'^a']y/I^ in the upper-right quadrant of the ms-m ^/2 
plane. The covariance is then given by. 



(5.19) 


and Pearson’s correlation coefficient is given by. 



(5.20) 


We note that the distributions of and 7773/2 can be found from equation (5.17) 
and are in fact not uniform, but peak at large values (e.g. the distribution for 7773/2 
is linear). Taking into account the full distribution would then lead to a stronger 
negative correlation, while including non-vanishing F-terms would effectively increase 
the radius of the permissible region in the 777 ^- 7773/2 plane, which would result in a 
smaller correlation. However, as we are not interested in the exact magnitude of 
this number but rather its sign, our simple discussion using uniform distributions 
suffices: we have shown that the continuous flux approximation generically gives 
r < 0, and for the special subset of points with 7773/2 -C [{2tiY a']y/Ll, it predicts 
r 0. Evidently, neither of these predictions explain the strong positive correlation 
we have observed in section 5.3, so we may now ask, why does the continuous flux 
approximation break down in the flux compactihcations that we consider? 

A possible explanation of the observed correlation would be to note that while the 
flux tadpole of Model 1, L* = 22, is larger than the number of flux cycles (2 x 4-1-2 = 
10), the hierarchy between these numbers is not necessarily large enough to justify the 
continuous flux approximation. Similarly, to hnd flux choices that satisfy the tadpole 
condition with some frequency, we have restricted the largest flux number to 5, which 
may well be much too small to justify the continuous flux approximation. For more 
general compactihcations than those considered in reference [47], signihcantly larger 
tadpoles can be found. We may then ask, does the relatively small hux tadpole drive 
the results of section 5.3? 

To address this question, we consider a hypothetical modihcation of the brane 
content of Model 1 that would give rise to a very large tadpole so that we may 
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(a) T = 5i, r = 0.92 


(b) T = lOi, r = 0.95 


Figure 7. The distribution of (m^, 771 , 3 / 2 ) at = —1.6 + 1.2i for 10,000 random flux choices 
drawn from a uniform distribution on [—50, 50] without imposing the tadpole condition. 

effectively take L* —)■ 00. Furthermore, we allow both positive and negative flux 
tadpoles and take the flux numbers to be uniformly distributed from —50 to +50. 
With these assumptions, the flux tadpole has no effect on the correlation coefficient, 
and the continuous flux approximation predicts r = 0. 

We show the resulting distribution of 7/13/2 and rris in Figure 7 for the reference 
point M* = —1.6 + 1.2* and r = 5i, 107 The resulting distribution is again strongly 
positively correlated with r = 0.92 and r = 0.95 - just as in the case plotted in 
Figure 6 for fluxes satisfying the tadpole condition. Thus, we conclude that the 
small tadpole of Model 1 does not affect or explain the correlation between rus and 
"13/2- 

The reason for the breakdown of the continuous flux approximation can instead 
be understood as follows. The large complex structure limit Im(M*) —)■ cx) is a so 
called ‘D-limit’ in which the vectors (11,11*, DiU, Djll*) no longer form a good basis, 
and Qfiux develops null directions [14, 17, 54]. Already in reference [14] it was argued 
that the continuous flux approximation should be expected to break down in the 
D-limit as the periods become comparable to the tadpole. 

In this paper we have seen that the continuous flux approximation breaks down 
quickly in the large complex structure expansion. This can be understood as a 
consequence of the fact that a single element (corresponding to the A°-cycle) quickly 
comes to dominate the norm of 11, and that the corresponding contribution to 77 ■ 11 
dominates the superpotential as a result. All scales of the resulting theory are then 
set by the fluxes on the cycle with all other fluxes only contributing by sub leading 
corrections. 

At weak string coupling, the contribution from the NS-NS fluxes to the superpo¬ 
tential are further enhanced by a factor of I/qs, and the RR fluxes give a subdomin- 
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ant contribution. A single real flux number then sets the scale of and m3/2, thus 
explaining the strong correlation in this case. 

We expect that the rapid dominance of the cubic terms in the superpotential 
is a generic feature of compactihcations with a large number of complex structure 
moduli. The number of cubic terms in the superpotential scales like as 

does the number of terms at quadratic order which come in sets of ~ terms, 

each proportional to one out of ~ flux choices. Thus, we expect that the results 
found in this paper should be applicable to compactihcations with a large number 
of complex structure moduli. 

5.5 Universality of random matrix theory and flux compactiflcations 

We are now ready to make two flnal observations regarding the applicability of ran¬ 
dom matrix theory to the statistics of the flux landscape. 

First, we note that the spectra that we have obtained analytically for M. and 
Ti in section 4 are not only simple and deterministic, but also have zero probability of 
appearing in the AZ-CI model, the ‘WWW’ model, or any straight-forward general¬ 
isation thereof. This can be understood as a direct consequence of the Vandermonde 
determinant that appears in the joint probability density of the eigenvalues. For 
example, the joint probability density of the N positive eigenvalues z/q of the AZ-CI 
ensemble is given by, 

/ 1 AT TV 

, . . . , Z/7v) OC YlWl- |exp ( In 

a<b \ a=l a=l 

where denotes the variance of the normally distributed independent matrix entries. 
The measure (5.21) clearly gives zero weight to configurations with degenerate ei¬ 
genvalues, such as the one we found in (4.22). 

While the Vandermonde determinant can be interpreted as giving rise to ‘ei¬ 
genvalue repulsion’, the observed spectra of flux compactiflcations considered in this 
paper are rather characterised by ‘eigenvalue attraction’ with peaked spectra and 
degenerate eigenvalues. We have shown that this spectacular difference arises due to 
the fast growth of a single element of the period vector, and is therefore particular 
to compactiflcations close to the large complex structure point. Such compactiflc¬ 
ations are interesting and much studied, but they only constitute a small fraction 
of all possible flux compactiflcations. Is it then the case that random matrix the¬ 
ory models are broadly applicable to flux compactiflcations on manifolds that are 
sufficiently distant from any D-limit? We here show that existing theorems in the 
Random Matrix Theory literature do not guarantee that the matrices A4 and H in 
the ‘flux landscape’ reach universal limits, but, we argue, there are good reasons to 
believe that the spectrum of Ad is well-described by the AZ-CI ensemble for generic 
compactiflcations with many moduli, as suggested in [16]. 
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The relevant RMT universality theorems can be understood as extensions of 
Wigner’s 1958 derivation of the limiting ‘semi-circle law’ for the eigenvalue density 
of an ensemble of symmetric matrices with independent and identically distributed 
entries [22], While that derivation relied on certain assumptions on the moments of 
the distribution of the independent matrix elements, it was independent of the details 
of the corresponding distribution and can therefore be regarded as ‘universal’. In the 
recent work [23], Wigner’s result was extended to matrices with correlated entries, 
and in [24] , similar theorems were derived for the remaining symmetry classes in the 
Altland-Zirnbauer classihcation, including the AZ-CI ensemble. Thus, reference [24] 
could provide the mathematical justihcation for asserting that the matrix M in flux 
compactihcations - despite string theory correlations - should roughly be described 
by the Cl matrix ensemble. 

References [23, 24] consider N x N matrices Mab in which some of the entries 
are in some way correlated, Mab ~ Mod- The exact nature of the correlation is 
unimportant, as long as the correlations are not too many and that the number 
of null vectors does not grow too fast with increasing matrix size. The three key 
assumptions are (in our notation) given by. 


Cl : 

max 

^ \ 

S {Mab ~ Mcd) | = 

KB,C,D J 

/ \ 

= 0{N^-^) , 

(5.22) 

C2 : 

max 

A,B,C 

f {Mab ~ McD)j = 

0 (W), 

(5.23) 

C3 : 

^ {Mab ~ Mbc) = 0{N^) 

A,B,C 

A^C 


(5.24) 


where e > 0 and. 


^ {Mab ~ Mcd) 


1 if Mab ~ Mqd , 
0 otherwise. 


(5.25) 


Assumption Cl ensures that the total number of correlations between all elements 
of a row and all elements of the matrix does not grow faster than 0 {N‘^~'^), while 
from assumption C2 it follows that no entry in the matrix has more than 0{N‘^) 
correlations with any full row of the matrix. Assumption C3 hnally states that 
the number of correlations between the R’th row and the R’th column (excluding 
elements related by the symmetries of Mab), when summed over B does not grow 
faster than 0 {N‘^). 

We now consider what these conditions mean for the matrix Ml in flux com¬ 
pactihcations that are described by the Gukov-Vafa-Witten superpotential (2.13) 
and the leading order Kahler potential (2.1).® We will assume that the independent 
components of Zab are Z^-i so that 6 {ZTi ~ Zrj) = 6{Z Ti 4) = sj. 

®We are very grateful to Kepa Sousa for discussions on this point. 
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The non-trivial correlations of M arise from the Hodge decomposition of the 
holomorphic three-form H and are given by eqnation (4.16) which we here may write 
as Zij = where C is an nnimportant constant. 

It will suffice for our purposes to consider the implications of Cl for A = r, even 
though this row does not necessarily have the maximum number of correlations of 
any row in the matrix. We then have, 

5 ~ Z,k) = Y, ^ ~ = Y ^ i^Ai) < 0{N^-^), (5.26) 

i,j,k ijik i,j,k 

for N = h]f + 1. Here 6 {xjki) = 0 if >Cjki = 0 and 5 {xjki) = 1 if >Cjki ^ 0. Condition 
C2 for C = T implies that, 


max 

*5 



r\j 


Zrk) 


= max 
hi 


= max 




< 0{N^). 


(5.27) 


Finally, condition C3 gives that. 


= + Z,Z, S{>Cjkp)S{>cikp) < 0{N‘^) . (5.28) 

a,h,c j,k j,l p,k 

a^c j^l 


Evidently, when interpreted in the context of flux compactifications the conditions 
(5.22)-(5.24) give rise to non-trivial conditions on the number of non-vanishing 
Yukawa couplings. 

For a set of compactification manifolds with an increasing number of complex 
structure moduli, the number of non-vanishing Yukawa couplings may in general 
scale with the total number of Yukawa couplings, i.e. ~ in violation of (5.26). 

Furthermore, the sum 5 {xijk) naturally scales linearly with h]f, in violation of 
(5.27) for e < 1. In addition, the last term of equation (5.28) may scale as fast as 
~ {h]YY, making also condition C3 inapplicable to flux compactifications. Thus, we 
see that the ‘universality theorems’ of [23, 24] do not automatically apply to flux 
compactifications. 

Obviously, this does not prove that random matrix theory is inapplicable to stud¬ 
ies of the flux landscape. It may be possible to construct a sequence of flux compacti¬ 
fications that in some region of the moduli space satisfy the conditions (5.26)~(5.28), 
thereby making the existing universality theorems applicable. Moreover, it may be 
that for correlations of the form (4.16), it is possible to derive stronger results than 
for the general correlations considered in [23, 24]. Indeed, by numerically simulating 
random matrices Ai that satisfy the string theory conditions (4.12) and (4.16) for a 
randomly chosen set of non-vanishing Yukawa couplings, we find that the eigenvalue 
density quickly becomes quite similar to the AZ-CI spectrum for large matrices. 
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We thus conclude that while existing universality arguments do not suffice to 
guarantee that the matrices M and "H in flux compactiflcations approach the Cl or 
‘WWW’ ensembles, respectively, it may well be possible to extend these arguments 
to the case that is particularly interesting in string theory. 

6 Conclusions 

We have considered four-dimensional effective supergravities arising in the low-energy 
limit of flux compactiflcations of type IIB string theory, and we have asked what are 
the spectra of the Hessian matrix "H and the matrix M that governs the critical 
point equation. By direct computation, we have found these spectra analytically in 
a subspace of the large complex structure limit in which the complex structure moduli 
all have the same phase. The resulting eigenvalue distributions are remarkably given 
by highly degenerate eigenvalues at integer multiples of |W| and 171^2 j independently 
of the details of the compactiflcation manifold or the flux choice. These results may 
thus be regarded as ‘universal’ for type IIB flux compactiflcations at large complex 
structure. 

By computing the spectra of "H and M numerically in explicit flux compactiflca¬ 
tions, we have found that the limiting spectra are quickly approached for Im(u*) > 1, 
and that over most of the sampled moduli space, 1713/2 is strongly linearly correlated 
with the scale of the supersymmetric moduli masses, m^. Such a correlation makes 
hierarchies of the form 7 / 73/2 irig highly infrequent in this region of moduli space. 

Our results imply that proposed random matrix theory models for "H and M are 
inapplicable at large complex structure, and we have furthermore argued that in 
more general type IIB flux compactiflcations, the number of correlations between 
matrix elements may grow too quickly for existing ‘universality theorems’ in the 
RMT literature to apply. It would be interesting to extend these theorems to include 
the particular relations that appear in flux compactiflcations. 

In this paper we have focussed on the moduli space dependence of eigenvalue 
distributions, and we have not specialised to the small subset of points that are 
supersymmetric vacua. It would be interesting to understand if the structure found in 
this paper can be found also among such vacua, or if on the contrary, supersymmetric 
vacua are preferably realised in the relatively small regions in which this structure is 
strongly broken. It would furthermore be interesting to extend the analysis of this 
paper to the full moduli sector of string compactiflcations. 

In sum, we have shown that in a particular region of the moduli space, the ‘flux 
landscape’ has much non-random structure that determines important aspects of the 
low-energy theory. Extending this study to broader regions of the moduli space is 
an important question for the future. 
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A Additional explicit examples of flux compactiflcations 

In this appendix, we present the results of numerical simulations of eigenvalue spectra 
in several additional examples of flux compactihcations. We begin with investigating 
the influence of the orientifold involution of the eigenvalue density in Model 1, finding 
that breaking the U2 O M3 symmetry by turning on general fluxes has no discernible 
effect on the eigenvalue distributions of M. and Ti. We will refer to this model as 
Model 1*. We then present the results for additional explicit flux compactihcations, 
namely Models 2-4 of [47] (here referred to by the same name) and the degree 18 
hypersurface in CP^ g g (here referred to as Model 5). 

The classical prepotentials are given in the large complex structure limit by 
[47, 48], 

F2 = —U^U2 — 3 MiM 2 ~ f "^2 ~ 2M1M2M3 — 2M2M3 — 2U1U2U4 — 2U2U4 — 2M2M3M4 + 2UiU2 
+M2M3 + U 2 U 4 + 1^2 + + |m2 + M3 + ^4 — *C(3)^ , (A.l) 

4 QQ Or^QO VO 070 

F3 = — gM]^ — 4MgM2 — 4 miM 2 — ^M2 — 4MgM3 — 8M1M2M3 — ^M2M3 — 4M1M3 — ^M2M3 

QQO QO QOTO 

—Mg — ^M 3 M 4 — 3M1M2M4 — IM2M4 — 3M1M3M4 — 3M2M3M4 — ^U^U4 — ^MiM4 
10 10 0 o o 

— 2M2M4 — 2U3U4 — M4M5 — 2M1M2M5 — M2M5 — 2M1M3M5 — 2M2M3M5 — M3M5 

—M1M4M5 — M2M4M5 — M3M4M5 + 2 ul + 4 miM 2 + 4miM3 + M1M4 + M2M4 + M3M4 


+ |m2 + |m2M3 + |m 3 + M1M4 + M2M4 + M3M4 + ^MiMg + ^U2U^ + 

+ yMi + ||m2 + 3 M 3 + |m 4 + Mg — , (A. 2 ) 

oQQioo o 72 

F4 = —u\ — ^U4U2 — ^MiM 2 — M4M3 — M1M2M3 — M4M4 — M1M2M4 — M1M3M4 — ^MgMg 

— 3M1M2M5 — ^MgMg — 2M1M3M5 — M2M3M5 — 2M1M4M5 — M2M4M5 — M3M4M5 — ^MiMg 

-|m 2M5 - M3M5 - M4M5 - ^Mg + |m4 + M1M2 + |miM3 + |miM 4 + ^MiMg + U2U5 

T^MsMg + |m5 + 3 Mi + |m 2 + M3 + M4 + §Mg - <( 3 )^ , (A. 3 ) 

^5 = -fw? - IM4M2 - |miM2 + |m4 + |miM 2 + ^Mi + |m2 - <( 3 )|f| . (A. 4 ) 
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The leading-order instanton corrections to the prepotentials are given by [47, 48],^ 


h 

h 

h 

h 


2g2i7™i 232e^*™2 2e^““^ -1-1886^*^“^ -t- 

^56e2m(«2+n3) + 5gg2*^(n2+n4) + _ _ _ ^ 

g2i7rtii _|_ 252e^*’^“® — 2e^*™^ -1- -|- e2*7r(Mi-i-'U2) _ 92526"^*™® 

-l-252g2“(“2+«3) _j_ g2i7r(ui-|-ii4) _|_ 2g2i7r(n4-|-ii5) 

252g2i7rMi I p2iTTU2 I 227ni3 , 2i7r4i4 , 2inus 


4 




(A.5) 

(A.6) 


_2g2*7r(“2+«4) _j_ 2ggg2i7r(Mi+n5) _|_ 

_j_35 •„2j7r?ii _ _Z_- 2i-KU2 _ 121_5 ,• 4i7r4ti i 45 4 ^447rU2 

27r3‘^ 87r3‘^ 167r3i67r3 

I n^ • 2i7r(iii+4t2) I 
“T 3 -r . .. , 


(A.7) 

(A.8) 


For Model 1* and Models 2-4 we consider a D3-tadpole contribution in the range 
0 < Qflux < 22 and uniformly sample the flux integers in the range [—5,5], while 
for Model 5 we use a D3-tadpole range of 0 < Qaux < 182 and a flux number range 
of [—25,25]. We here present the numerical results for the spectra of M. and H 
computed in the same cases as was done for Model 1 in Figures 2 and 3. 
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